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The excited states of a charged particle interacting with the quantized electromagnetic field 
and an external potential all decay, but such a particle should have a true ground state — one 
O |. that minimizes the energy and satisfies the Schrodinger equation. We prove quite generally that 

fH . this state exists for all values of the fine-structure constant and ultraviolet cutoff. We also show 

, the same thing for a many-particle system under physically natural conditions. 

a, 

: 1 INTRODUCTION 

X 

An established picture of an atom or molecule is that even in the presence of a quantized radiation 
field there is a ground state. The excited states that exist in the absence of coupling to the field 
are expected to melt into resonances, which means that they eventually decay with time into the 
ground state plus free photons. This picture has been established by Bach, Frohlich and Sigal in 
H for sufficiently small values of the various parameters that define the theory. Here we show 
that a ground state exists for all values of the parameters (including a variable ^-factor) in the one 
particle case and, under a physically appropriate assumption, in the many-particle case. 

We know that the Hamiltonian for the system is bounded below, so a ground state energy 
always exists in the sense of being the infimum of the spectrum, but the existence of a genuine 
normalizable solution to the eigenvalue equation is a more delicate matter that has received a great 
deal of attention, especially in recent years. A physically simple example where no ground state 
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exists (as far as we believe now) is the free particle (i.e., particle plus field). In the presence of an 
external potential, however, like the Coulomb potential of a nucleus, a ground state should exist. 

The difficulty in establishing this ground state comes from the fact that the bottom of the 
spectrum lies in the continuum (i.e., essential spectrum), not below it, as is the case for the usual 
Schrodinger equation. We denote the bottom of spectrum of the free-particle Hamiltonian for 
N particles with appropriate statistics by E°(N). The 'free-particle' Hamiltonian includes the 
interparticle interaction (e.g., the Coulomb repulsion of electrons) but it does not include the 
interaction with a fixed external potential, e.g., the interaction with nuclei. When the latter is 
included we denote the bottom of the spectrum by E v (N). It is not hard to see in many cases that 
E v (N) < E°(N), but despite this inequality E v (N) is, nevertheless, the bottom of the essential 
spectrum. The reason is that we can always add arbitrarily many, arbitrarily 'soft' photons that 
add arbitrarily little energy. It is the soft photon problem that is our primary concern here. 

The main point of this paper is to show how to overcome this infrared problem and to show, 
quite generally for a one-particle system, that a ground state exists for all values of the particle 
mass, the coupling to the field (fine-structure constant a = e 2 /hc), the magnetic g-factor, and 
the ultraviolet cutoff A of the electromagnetic field frequencies, provided a bound state exists 
when the field is turned off. This result implies, in particular, that for a fixed ultraviolet cutoff 
renormalization of the various physical quantities will not affect the existence of a ground state. Of 
course, nothing can be said about the limit as the cutoff tends to infinity. We also include a large 
class of interactions much more general than the usual Coulomb interaction. 

The model we discuss has been used quite frequently in field theory. In its classical version it 
was investigated by Kramers [O] who seems to have been the first to point out the possibility of 
renormalization. The quantized version was investigated by Pauli and Fierz [p4|| in connection with 
scattering theory. Most importantly, it was used by Bethe M to obtain a suprisingly good value 
for the Lamb shift. 

Various restricted versions of the problem have been attacked successfully. In the early seventies 
Frohlich investigated the infrared problem in translation invariant models of scalar electrons coupled 



to scalar bosons [[b|]. In [14] he proved that for an electron coupled to a massive field a unique 
ground state exists for fixed total momentum. External potentials were not considered in these 
papers. 

The first rigorous result on the bound state problem, to our knowledge, is due to Arai and 
concerns one particle confined by an x 2 -potential, the interaction with the photons being subject 
to the dipole-approximation. For this model, which is explicitly solvable, Arai proved existence and 
uniqueness of the ground state ||. Later, Spohn showed by perturbation theory that this result 
extends to perturbed x 2 - potentials [26]. No bounds on the parameters were needed to obtain these 
results but the methods oviously do not admit extensions to more realistic models. 

Bach , Frohlich and Sigal 0, |?]] initiated the study of the full nonrelativistic QED model (the 
same model as considered in the present paper) under various simplifying assumptions. In || the 
existence of a ground state in this model for particles subject to an external binding potential was 
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proved for oA small enough. The main achievement of this paper is the elimination of the earlier 
simplifying assumptions, especially the infrared regularization. This is the first, and up to now the 
only paper where a 'first principles' QED model was successfully analyzed, but with a restricted 
parameter range. A weaker result for the same system but with simplifications such as infrared 
regularization were independently obtained by Hiroshima by entirely different methods; he also 
showed uniqueness of the ground state, assuming its existence [17]. 

In a parallel development Arai and Hirokawa ||, f|], Spohn |2^], and Gerard [|l^] investigated the 
ground-state problem for systems similar to the one of Bach et al in various degrees of mathematical 
generality. Arai and Hirokawa proved existence in what they called a 'generalized spin-boson model'. 
If specialized to the case of a non-relativistic A r -particle system interacting with bosons, their result 
proves the existence of a ground state for confined particles and small a. This result was extended 
in a recent preprint @ to account for non-confined particles and systems with the true infrared 
singularity of QED. Concrete results in the infrared-singular case concern only special models, 
however, such as the Wigner-Weisskopf Hamiltonian which describes a two-level atom. Hirokawa 



continued this work in a recent preprint [16]. Spohn and Gerard both proved existence of a ground 
state for a confined particle and arbitrary coupling constant, the result of Gerard being somewhat 



more general [28, 15 1. 

For the existence of a ground state in the case of massive bosons, which is a typical intermediate 
result in the cited works, a short and elegant proof was given by Derezinski and Gerard [11, 15] for 
the case of linear coupling, in which the A 2 term is omitted, and a confining potential. Some of 
their ideas are used in our paper. 

The Hamiltonian for A^-particles has four parts which are described precisely in the next section, 

H v = T + V + I + H f . (1) 

The dependence of H v on N is not noted explicitly. The first term, T = Ylj=i ^ s the kinetic 
energy with 'minimal coupling' in the Coulomb gauge (i.e., p is replaced by p + y/a A, where A is 
the magnetic vector potential satisfying divyl = 0), V(X), with X = (xi,X2, ...,xn) is the external 
potential, typically a Coulomb attraction to one or more nuclei. In any case, we assume that V is 
a sum of one-body potentials, i.e., 

N 

V{X) = Y,v{x 5 ) . 
j=i 

The particle-particle interaction, /, has the important feature that it is translation invariant 
and, of course, symmetric in the particle labels. Both / and V could be spin dependent, but we 
shall not burden the notation with this latter possibility. Typically / is a Coulomb repulsion, but 
we do not have to assume that / is merely a sum of two-body potentials. The only requirements 
are: (1) the negative part of v vanishes at infinity; (2) the negative part of I satisfies clustering, i.e., 
the intercluster part of J_ tends to when the spacing between any two clusters tends to infinity; 
(3) V- + 2_ are dominated by the kinetic energy as in (||). Another assumption we make (in the 
A^-particle case) is that there is binding, as described below. 
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The natural choice for the kinetic energy is the 'Pauli operator' T = (jp + \fa A(x)) 2 + \fa a ■ B, 
but we can generalize this to include the case of the usual kinetic energy T = (p + \faA(x)) 2 by 
introducing a 'g-factor', g £ R. Thus, we take 

T= (p + V^0z0) 2 + fVa^-£ . (2) 

Note that T is a positive operator if < g < 2 and may not be otherwise. Nevertheless, the 
Hamiltonian is always bounded below because of the ultraviolet cutoff we shall impose on the A 
field, which implies that (g/2)^/aa ■ B + Hf is always bounded below. 

We believe that the 'relativistic' operator T = \p + ^/aA(x)\, presents no real difficulty either, 
but we do not want to overburden this paper with a lengthy proof. This problem is currently under 
investigation. 

A model that is frequently discussed is the 'Pauli-Fierz' model, but it is not entirely clear how 
this is defined since several variants appear in |24j]. One version uses T = (p + ^/aA(x)) 2 + yfaa ■ B, 
which is one of the models under consideration. Another variant uses a linearized version of this 
operator, T = p 2 + 2^fap ■ A(x) + ^fa a • B or T = p 2 + 2yfap • A(x). These variants are not gauge 
invariant and, therefore, depend on the choice of gauge for A. Our method is applicable to these 
linearized models in some gauges, but not in others. We omit further discussion of this point since 
these variants are not the most relevant ones for quantum electrodynamics. 

There is one important point that as far as we know, has not been mentioned in the QED 
context. This is the binding condition. Our proof of the existence of a ground state uses, as input, 
the assumption that 

E V (N) < E V (N') +E°(N- N') for all N' < N. (3) 

Our work can be generalized (but we shall not do so here) to the case in which the external 
potential is that of attractive nuclei and these nuclei are also dynamical particles. Then, of course, 
it is necessary to work in the center of mass system and then (||) must be replaced by the condition 
that E v (N) is less than the lowest two-cluster threshold. While this condition, or (g) in the static 
case, are physically necessary for the existence of a ground state, the validity of these conditions 
cannot be taken for granted. 

We prove inequality (|3|) for one particle (N = 1) quite generally, using only the assumption that 
the ordinary Schrodinger operator p 2 +v has a negative energy ground state. This certainly holds for 
the Coulomb potential. Indeed, one could expect, on physical grounds, that there could be binding 
even if p 2 + v is has no negative energy bound state, because the interaction with the field increases 
the effective mass of the particle — and hence the binding energy. The same argument shows that 
when there are iV particles at least one them is necessarily bound, i.e., E v (N) < E°(N). When we 
consider more than one particle, we are not able to show (^) for all N 1 > 1, even if ^ p 2 + / + V 
has a ground state. 

In the Coulomb case, it is possible to show (but we shall not do so here) that condition (|3|) is 
satisfied if the nuclear charge Z is large enough. The basic idea is that if breakup into two groups, 
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one of them with N' particles close to the nucleus and the second consisting of N — N' particles 
far away occurs, then there will be an attractive Coulomb tail acting on the separated particles 
at a distance R away with net attractive potential (Z — N')/R. However, to localize one of these 
particles within a distance R of the nucleus will require a field energy localization error of the order 
of C/R, by dimensional analysis arguments. If (Z — N') > C then the energy can be decreased by 
bringing one of the unbound particles close to the nucleus. 

Section || introduces the precise definition of our problem and the main result Theorem ^jj . 

In Section 3, we show how to prove that E v '(1) < E°(l). More generally, E V (N) < E°(N) if 
v < 0. 

Our strategy to establish a ground state is the usual one of showing that a minimizing sequence 
of trial vectors for the energy actually has a weak limit that, in fact, is a minimizer. The problem 
here is that one can easily construct minimizing sequences that converge weakly to zero by choosing 
vectors with too many soft photons. To avoid this we take a special sequence. 

To define this sequence we first consider an artificial model in which the photons have a mass, 
i.e., oj(k) = \Jk 2 + m 2 . Here there is no soft photon problem and we show in Section 4 that this 
model has a ground state $ m . 

In Section 5 we show that as m — > the <& m sequence is minimizing. Then in Section 6 we use 
the Schrodinger equation for $ m to deduce certain properties of & m which we call infrared bounds. 
One of these was proved in ||] but we need one more, which is new. 

With these bounds we can show in Section 7 that $ m has a strong limit as m — > 0, which is a 
minimizer for H v . 

Acknowledgment: We thank Professor Fumio Hiroshima for a useful correspondence concerning 
equation (|53|) and for sending us his preprint [^] . 

2 DEFINITIONS AND MAIN THEOREM 

The Hamiltonian for N particles interacting with the quantized radiation field and with a given 
external potential V(X), with X = (xi,X2, ...,xn) and xj S M 3 , is 

N 

RV = Y,{ & + ^A( Xj )) 2 + • B( Xj )] + V(X) + I{X) + H f . (4) 

The unit of energy is Mc 2 /2, where M is the particle mass, the unit of length is £ c = 2h/Mc, 
twice the Compton wavelength of the particle, and a = e 2 /he is the dimensionless "fine structure 
constant" (=1/137 in nature). The electric charge of the particle is e. The unit of time is the time 
it takes a light wave to travel a Compton wavelength, i. e., the speed of light is c = 1. 

The operator p = — iV, while A is the (ultraviolet cutoff) magnetic vector potential (we use the 
Coulomb, or radiation gauge). The unit of A 2 {x) is Mc 2 /2£ c . The magnetic field is B = cur LA and 
the unit of B is o?l 2 times the quantity M 2 e 3 c/4ft? , which is the value of B for which the magnetic 
length (hc/eB) 1 / 2 equals the Bohr radius 2h 2 /Me 2 . 
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The reader might wonder why we use these units, which seem to be more appropriate for a 
relativistic theory than for the nonrelativistic theory we are considering. Why not use the Bohr 
radius as the unit of length, for example? Our reason is that we want to isolate the electric 
charge, which is the quantity that defines the interaction of matter with the electromagnetic field, 
in precisely one place, namely a. 'Atomic units' have the charge built into the length, etc. and we 
find this difficult to disentangle. 

The Hilbert space is an appropriate subspace of 

where T is the Fock space for the photon field. We have in mind Fermi statistics (the antisymmetric 
subspace of ® N L 2 (M. 3 ; C 2 )) and the C 2 is to accomodate the electron spin. We can also deal with 
"Boltzmann statistics" , in which case we would set g = and use ® N L 2 (M. 3 ), or bose statistics, in 
which case we would set g = and use the symmetric subspace of (g) L 2 (R 3 ). These generalizations 
are mathematically trivial and we do not discuss them further. 

For our purposes we assume that for every e > there exists a constant a(e) such that the 
negative part of the potentials, V—(X) and I-(X), satisfy 

N 

V. + I-Ke^pj + aie) (5) 



as quadratic forms on H. 
The vector potential is 



A{x)=Y J I —[ex(k)a x (k)+e x (-k)al(-k)]e ik - x d 3 k (6) 



1 

A=l,2-l fc l< A VT^I 

where the operators a x ,a* x satisfy the usual commutation relations 

[a x (k),al(q)} = 5(k-q)5 Xl v , [a x (k), a u (q)] = 0, etc (7) 

and the vectors £\(k) are the two possible orthonormal polarization vectors perpendicular to k. 
They are chosen for convenience in ( |59| , |60|) . 

The number A is the ultraviolet cutoff on the wavenumbers k. Our results hold for all finite 
A. The details of the cutoff in (^) are quite unimportant, except for the requirement that rotation 
symmetry in A;-space is maintained. E.g., a gaussian cutoff can be used instead of our sharp cutoff. 
We avoid unnecessary generalisations. 

The field energy Hf, sometimes called dT(uj) is given by 

H f =J2 [ co(k)a* x (k)a x (k)d 3 k (8) 

A=l,2 

The energy of a photon is uj(k) and the physical value of interest to us is 

u(k) = \k\, (9) 
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in our units. Indeed, any continuous function that is bounded below by const. \k\ for small \k\ is 
acceptable. In the process of proving the existence of a ground state for H we will first study the 
unphysical 'massive photon' case, in which 

u m (k) = \jk 2 + m 2 (10) 

for some m > 0, called the 'photon mass'. 

In the remainder of this paper, unless otherwise stated, we shall always assume that there is no 
restriction on a, A and g and that uj{k) can be either as in (|9|) or as in fllpp. 

By Lemma A. 5 we see easily that H v is bounded below for all values of the parameters, including 
m = 0. Thus, H v defines a closable quadratic form and hence a selfadjoint operator, the Friedrich's 
extension. We denote this extension again by H v . 

Our main theorem is 

2.1. THEOREM (Existence of a ground state). Assume that the binding condition and 
the condition (j^) hold. Then there is a vector <J> in the N -particle Hilbert space 7i such that 

H v $ = E V (N)<S> . (11) 



3 UPPER BOUND 

We shall prove the binding condition (|3|) for one particle and, with an additional assumption, for 
the N particle case as well. This is that if A-particles are present then at least one of them binds. 

As we mentioned before, at least in the single particle case, that our requirement that the 
system without the radiation field has a bound state is somewhat unnaturally restrictive, since one 
expects that the radiation field enhances binding; this has been shown to be true in the 'dipole', 
or Kramers approximation |19| . We are able to show in the one-particle case, that the photon field 
cannot decrease the binding energy. It is quite possible that there could be binding even when 
the operator p 2 + v does not have a negative energy state, but we cannot shed any light on that 
question. 

For the one-particle case the situation is less delicate than the A^-particle case. 

3.1. THEOREM (Binding of at least one particle). Assume that the one-particle Hamilto- 
nian p 2 + v(x) has a negative energy bound state with eigenfunction (f)(x) and energy — eo- Then, 

E v (l)<E°(l)-e , (12) 

i.e., binding continues to exist when the field is turned on. 

For the N -particle case we make the additional assumption that v(x) < for all x. Then, 

E v (N)<E°(N)-e , (13) 

i.e., at least one particle is bound. 
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Proof. It suffices to prove that E V (N) < E°(N)+e- e for all e > 0. There is a normalized vector 
F € TL such that (F,H°F^J < E°(N) + e. (F is antisymmetric according to the Pauli principle.) 
We use the notation (•, •) to denote the inner product in Fock space and spin space. Then we can 
write (F, H°F) = J G(X)d m X with 

G(X) = 

N 

^{((-iVj + V^A( Xj ))F, HV, + JEA( Xj ))F)(X) + ^ (s/2) <F, ^ ■ B( Xj )F)(X)} 
3=1 

+ (F,(I + H f )F)(X) . (14) 



As a (unnormalized) variational trial vector we take the vector ip = YljLi ^ > { x jY 



1/2 



F. Recall 



that 4>(x) > since <f> is the ground state of p 2 + v. We also recall the Schwarz inequality 



Ef=i<K^)v<^) 



2_,=i 



1/2 



N 

3=1 



(15) 



Using (|l5|), integration by parts, and the fact that (ft satisfies the Schrodinger equation (y? + 
v)<p = —eo4>, we easily find that 



(V, [H v -(E°(N)+e-e )] ^) 

< J {G(X) - (E°(N) +s)(F,F)(X)}Y,<P(x j ) 2 d 3N X 
+ jY,v(^)Hxj) 2 (F,F)(X)d 3N X . 



N 



3=1 



(16) 



When N = 1 the last term in ([lq ) is not present so no assumption about the potential v is 
needed. When iV > 1 we can omit the last term because it is negative by assumption. 

Now, by the M 3 -translation invariance of Hq, for every y € M 3 there is a 'translated' vector F y 
so that G{X) -» G(X + (y, ...,y)) and (F y ,F y )(X) = (F,F)(X + (y, ...,y)). (This is accomplished 
by the unitary operator on Ti. that takes Xj Xj + y for every j and a x (k) — > exp(ik ■ y)a x (k).) 



Thus, if we denote the quantity in {} in (16) by W(X), and if we define ip y by replacing F by F y 
in the definition of ip, we have 



„ N N 

< I W(X + (y 1 ...,y))^<P(^) 2 d 3N X= I W(X) ]T <f>( Xj - y?d 3N X . 

3=1 3=1 



(17) 



Note that / Sl(y)dy <N J W(X)d 3N X. But / W(X)d 3N X = (F, (H° - E°(N) - e) F) and this 
is strictly negative by assumption. Hence, for some y £ M 3 we have that £l(y) < and thus tp y ^ 0, 
which proves the theorem. □ 



Remark: [Alternative theorem] 
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It may be useful to note, briefly, a different proof of Theorem 3.1, for long range potentials v(x), 
such as the attractive Coulomb potential —Z/\x\, which shows that the bottom of the spectrum of 
H lies strictly below E°. Unfortunately, this proof does not show that the difference is at least 
eo- We sketch it for the one-body case. Using the notation of the proof above, the first step is 
to replace F by F R = u{x\/R)F where u is a smooth function with support in a ball of radius 1. 
One easily finds that (F R , H F R ) / (F R , F R ) = E° + e + c/R 2 , where c is a constant that depends 
only on u and not on e and R. On the other hand (F R , V F R ) / (F R , F R ) < —Z/R, to use the 
Coulomb potential as an example. To complete the argument, choose R = 2c/ Z and then choose 
e = c/R 2 . What we have used here is the fact that localization 'costs' a kinetic energy R- 2 , while 
the potential energy falls off slower than this, e.g., R . 



4 GROUND STATE WITH MASSIVE PHOTONS 

As we emphasized in the introduction, not every minimizing sequence converges to the minimizer 
for our m = problem, i.e., with uo(k) = \k\. The situation is much easier for the massive case 
(|l0|). The Hamiltonian in this case is given by (||) and Hf is given by (||) with (|l0|). To emphasize 
the dependence on m we denote this Hamiltonian and field energy by and -ffy(m), respectively. 
Likewise, E v (m,N) and E°(m,N) denote the mass dependent energies, as defined before. 

We emphasize that the vector potential is still given by (|6|), but we could, if we wished, easily 
replace \k\^l 2 in (§) by (k 2 + m 2 ) -1 / 4 . 

It will be shown in this section that has a ground state. More precisely we prove 

4.1. THEOREM (Existence of ground state). Assume that for some fixed value of the ul- 
traviolet cutoff A there is binding for the Hamiltonian H^, i-e., E v (m,N) < S (m, N) where 
£ y (m, iV) = min{E v (m,N') +E°(m,N - N') : all N' < N} is the 'lowest two-cluster threshold'. 
Then E(m, N) is an eigenvalue, i.e., there exists a state *&rn i^ such that Fl^^m — E(rn, A^)$ m . 

Proof. Let us first show that it suffices to prove that for any normalized sequence ^ J , j = 1, 2, 
(not necessarily minimizing) tending weakly to zero 

liminf(^',i/ y ^') > E v (m,N) . (18) 

j—*oo 

To prove this let $ J be some minimizing sequence, i.e., assume that 

11^11 = 1, (19) 

and that 

(&,HV&)^E v (m,N) . (20) 

By the Banach Alaoglu Theorem we can assume that this sequence, as well as the sequence 
converge weakly in the sense that for any ^ € TL with (^,H^ l 1 ^) < oo we have that 

->(*,iZ£$ m ) , (21) 
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where <E> m is the weak limit of f» J . Our goal is to show that ($ m ,iJ^$ m ) = E (m,N) and that 



$m = I- 



Write ® J = $ m + VF. Obviously \I>- ? as well as il^*- 7 go weakly to zero. Thus 

= lim ($ j , (F y - E y (m, iV))^') 
= lim (($ m + {Hi - E v (m,N))(<Z> m + &)) 

j-*oo 

= lim (H% - E v (m, N))&) + (* m , (H v m - E v (m, N))$ m ) 

where we used that the cross terms vanish. Since Hl-E v (m, N)>0 this shows that $ m minimizes 
the energy, and, furthermore, that 

> lim (&, (Hi - E v (m,N))&) > <51iminf \\&\\ 2 

j^oo j—*oo 

for some positive constant 5. The second inequality is trivial if liminfj^oo H^'H 2 = and otherwise 
follows from our assumption (]i~8|). This proves that ^ converges strongly to zero along a subse- 
quence, which implies that [|$rn|| = 1- Hence <3? m is a normalized ground state. Thus, it suffices to 



prove (18) 



The steps that lead to a proof of (|18| ) are quite standard. The only difficulty is that one has to 
localize in Fock space, which we describe first. We follow [11] with some necessary modifications 
and some simplifications. 

Recall that, when the a x operators are viewed in x-space 

a x (f) a\(g) : T ^ T , (22) 

they obey the commutation relations 

[a x (f),a* x (g)) = [ J(x)g(x)d 3 x =: (f,g) . (23) 



Consider now two smooth localization functions j\ and j% that satisfy j\ + j\ = 1 and j\ is 
supported in a ball of radius P. The first derivatives of j% and j2 are of order 1/P. 
The operators 

cx(f) = a x (j 1 f)®l + l®a x (j 2 f) ,c* x (g)=a* x (j 1 g)®I + X®al(j 3 g) (24) 

act both on the space T ® T . Note that 

[cx(f),c* x (g)] = (f,g) . (25) 

Thus, these new creation and anihilation operators create another Fock space T l that is a subspace 
of T ® T and is isomorphic to the old Fock space T . Hence, there exists a map 

U :T (26) 
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that is an invertible isometry between Fock spaces. It is uniquely specified by the properties 

a* = U*c*U , (27) 

and the vacuum in T is mapped to the vacuum in T % T . 

The map U* is defined on T l only, but we can extend it to all of T ® T by setting U*F = 
whenever F 6 T <S> T is perpendicular to T l . In other words £7* is a partial isometry between 
Fockspaces where U*U = I on J 7 , and where [/[/* is the orthogonal projection onto T . We 
continue to denote the extended map by U*. 

Let 4> an d <f> be smooth nonnegative functions, with (j) 2 + cj) =1,0 identically one on the unit 
ball, and vanishing outside the ball of radius 2. Set (p R (X) = <j>(X/R). It is a standard calculation 
to show that for any ^ with finite energy 

(% H^S) = (<f> R H>, Hl^) + ffS$) - (*, (V0 fl ) 2 *) - (tf , (V^) 2 *) . (28) 



The last two terms in (|28| ) are bounded by const./ R 2 . 

One goal will be to show that for any ^ with finite energy 

0*17* ® X + J ® #/} ) + o(l) • (29) 

The error term o(l) vanishes as both R and P go to infinity and depends otherwise only on the 
energy of Notice that the invertible map U depends on the cutoff parameter P as well. (|29|) will 



be proved in Lemma A.l in the Appendix. The intuition behind the estimate ( |29| ) is that localized 
electrons interact only weakly with far away photons. Those photons are described solely by their 
own field energy. 



An immediate consequence of ( p9[ ) is the estimate 

(%<(> R H^<I> R 9) > (E v (m,N)+m)\\<p R n 2 -rn(<p R y,U*l®P2U0 R y) + o(l) , (30) 

which is obtained by noting that the field energy in the second factor can be estimated from below 
by 

H f >ml- mP 2 , (31) 
where P 2 is the projection onto the vacuum of the second factor of J- ® T . 



In a further step we prove in Lemma A. 3 that the sequence 



(fa**, U*I ® IW(f>R* j ) (32) 



as j — > oo. 



Returning to (|28|), using Corollary A. 2 we have that 



■^,^^)>i; v (m,JV)||^|| 2 -o(l) , (33) 
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with o(l) going to zero as R — > 00. Roughly speaking (ft forces some of the particles to be far 
away from the origin. Any such particle configuration can be described by two clusters with no 
interaction between them. In particular the interaction between these clusters via the radiation 
field is turned off. This means that each cluster carries its own field energy. To prove this the 
localization in Fock space is used. Moreover the cluster that is far away from the origin does not 
interact with the external potential although the repulsion among its particles is still present. 
To summarize, by combining fl30|) , (|32|) and (]33|) we have proved that 

liminf^',^* 3 ') > (E v (m, N) + S) + o(l) (34) 

where 

5 = min{m, S v (m, N) - E v (m, N)} , (35) 

and o(l) tends to zero as R — ► 00 and P — > 00. □ 



5 A MINIMIZING SEQUENCE 

We consider the Hamiltonian defined in (||) with field energy Hf(m) defined using u) m {k) = 
\/k 2 + m 2 . Our main goal here is Theorem 5.3, which shows that the ground states of the m > 



problem form a minimizing sequence for the m = problem. 

5.1. THEOREM (E v (m,N) converges to E v (0,N)). As m -> 0, 

E v (m,N) -» E v {0,N) and E°(m, N) -> E°(0, N) (36) 

Proof. First, note that > HY., > HY if m > m! > 0, because u> m has this same monotonicity 
property. Therefore, for any sequence of m — > 0, E v (m, N) is monotonically decreasing and has a 
sequence-independent, finite limit, which we call E* , and we note that E* > E v (0,N). To prove 
the opposite, namely E* < E v (0, N), we shall prove that E* < E v (0, N) + 2e for every e > 0. 

Let $ G W be normalized and such that ($, $) < £^(0, iV) + e. 

We note that < Hq + mAf, where Af is the number operator 



Af= [ a* x (k)a x (k)d 3 k. (37) 

A=l,2 ^l fc l <A 

Thus, if we use $ as a variational function for we have E v (m,N) < E v {0,N) + e + 
m ( C &,AA $), and our goal is accomplished provided that $ can be chosen so that ($, N <3?) < cxd, 
in addition to ($, -ffg 7 < E v (0, iV) + e. If a way can be found to modify $ to another vector <5 
so that (<fr,N &j < 00, in addition to ($>,H^ < E v (0,N) + 2e the proof will be complete. 

A suitable choice is $ = II„<I>, with Tl n being the projector onto the subspace of 7i, with n or 
fewer photons, i.e., 7i n = ® Ar L 2 (IR 3 ; C 2 ) ®T< n , for an appropriately large n. It is easy to see, with 
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the help of Lemma [AJ that (II n $, [V + I]n„$) -► ($, [V + /]$), that (n n $, #/]!„$) -> ($, #/$), 
and (Iln^jplnn^) — > (<£,p|<J>) as n — > oo. The following Lemma |5.2| , shows that the other terms 
converge as well. The same proof works for E°(m, N). □ 

5.2. LEMMA (Finite photon number approximation). Let $ G H be such that (&,Hf$) < 
oo. Denote by Tl n the projection in 7i onto states with photon number less than or equal to n, i.e., 
n n is the projection onto the subspace 7i n = ® Ar L 2 (]R 3 ; C 2 ) (g) J-< n C TL. Let <J> n = n n <3?. Then we 
have the following strong convergence as n — ► oo (in addition to <3? ra — > <3? ) 

Bfo)*,, -> Bfo)* . (38) 

Proof. Write A = D + D* where D contains the annihilation operators and D* contains the creation 
operators. We omit (xj) and we omit the vector index of A for simplicity. Since Hf&) < oo we 
learn from Lemma |A.4j that and D*§ are in 7i. Since IT n _i — ► I strongly, 

The same holds if -D is replaced by D* and n — 1 by n + 1. This proves the statement for A. The 
statement for B is proved in the same way. □ 



The following are two corollaries of Theorem 5T| . The first is the input for Section |g. The 



second is important for showing that it is only necessary to state the 'binding' condition (0) once 



(for m = 0). It is a trivial consequence of Theorem 5.1. 



5.3. THEOREM (Minimizing sequence). Suppose that mi > mi > ... > is a sequence 
tending to zero and suppose that <&j, for j = 1,2,... is an approximate minimizer for H^. in the 
sense that 



Sj = Hi. - E v (mj, N) - as j - oo . (39) 
Then $i, <&2> ••• * s a minimizing sequence for Hq . 

Proof. E v ( mj , N) + Sj = (*,-, Hi >(*,-, i# $,) > E v (0, N). 



□ 

5.4. LEMMA (Binding without mass implies uniform binding with mass). Assume ||) ; 
i.e., assume that 

E v (0, N) < mm{E v (0, N') + E°(0, N — N') : N' < N} — 2e . 
Then, for all sufficiently small m, 

E v (m, N) < min{E v (m, N') + E°(m, N — N r ) : N' < N} — e . 



Proof. Each of the various energies converges as m — > by Theorem 5.1. □ 
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6 TWO INFRARED BOUNDS 

We have seen that any sequence of (approximate) minimizers for the problem (m > 0) is 
a minimizing sequence for the HY problem. A crucial point was the possibility of finding an 
approximate minimizer for the HY problem that has a finite expectation value for the total photon 
number operator J\f . 



We also know by Corollary 5.4 and Theorem 4.1 that the HjL problem has a ground state, $ 



and in this section we shall prove two theorems about the soft photon behavior of $ m . The first, 
Theorem 6.1, is about the photon number a x (k)a x (k) and is based on a method of ||. The second, 



Theorem |T^, which has no antecedent we are aware of, is about the derivative of a\ (k)a x (k) with 
respect to k. 

6.1. THEOREM (Photon number bound). Assume that there is binding, i.e., E v (m,N) < 
T, v (m, N). Assume that $ m is a normalized ground state for the many-body Hamiltonian H^, 
m > 0. Then 

Pa 

($ m ,a* x (k)a x (k)$ m ) < —-(l+g 2 ) X A(k), 



where P is a finite constant independent of& m , g, a, and depends on m only via the binding energy 
T, v (m, N) — E v (m, N) > and of course on A . The function x\(k) is the characteristic function 
of the ball of radius A. 



Remark: We have proved in Corollary 5.4 that the binding energy is a uniformly positive 



function of m for small m, if the binding energy at m = is not zero. Therefore, Theorem 6.1 
implies that the number operator M in (|37|) is uniformly bounded for small m. 



The proof of Theorem 3.1 will be based on the following lemma about exponential decay of 



eigenf unctions, concerning which there is a vast literature (see HH 0, |2§ §). We do not strive 
at all to get the best exponential decay constants. For us the only really relevant goal is a decay 
estimate that depends only on T, v (m, N) — E v (m, N) > 0, but not otherwise on m. 

6.2. LEMMA (Exponential decay). Let be the N-body Hamiltonian in (Qj and let Q m be 

a groundstate wave function, which necessarily satisfies the Schrddinger equation 

H^ m = E v (m,N)^> m . (40) 

We assume that T, v (m,N) - E v (m,N) > and choose (5 > with (3 2 < T, v (m,N) - E v (m,N). 
Then 

ll«*K/TOM' < c(i + ^^J^^V II-M 2 (4i) 

where the constant C does not depend on m. 



The strategy of the following proof is probably due to Agmon Q . We learned it from 1 2C ] 
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Proof. Let G(X) any smooth, bounded function on M. 3N with bounded first derivative. We easily 
compute 

[ [Hl-E v (m,N), G], G]=-2\VG\ 2 . (42) 
We use the Schrodinger equation ( |40| ) to compute 



(G$ m , [Hi - E v (m,N)}G$ m ) = -~ ($ m , [ [H% - E(m), G] , G] * m ) = ($ m , |VG| 2 $ m ) . 

(43) 

Now we choose G to be 



G(X) = X (X/R)exp[f(X)}, where 

nx) r w 



(44) 



,l + e|X| 

and where < x ^ 1 is a smooth cutoff function that is identically equal to 1 outside the ball of 
radius 2, and identically zero inside the ball of radius 1. We let e — > at the end. 
Next, we calculate 

|VG| 2 = |V X | 2 e 2/ + 2V X • Vfe f G + |V/| 2 G 2 , 

and note that the first and second terms are compactly supported in M. 3N and each is bounded by 
a constant C that depends on (3 and R. 

Returning to (^), we obtain, after rearranging terms, 

(G<& m , « - £ y (m, AT) - |V/| 2 )) G$ m ) < C\\1> m \\ 2 . (45) 



Since |V/| < (3 we know by Corollary A. 2 that 



X«-^(m,iV)-|V/| 2 )x 
> (S y (m,AT) -E v (m,N) - (3 2 - o(l)) X 2 (R oo) . (46) 



In conjunction with (]45|) this shows that 



l^d, ||2 < ^ 11$ 

1 m " " S^(m, AT) - #(m, AT) - /3 2 11 r 



for i? large enough. After letting e — ► by monotone convergence a similar bound with G replaced 
by xexp(/?|X|) is obtained. □ 



Proof of Theorem 6A_. This proof is a slight modification of the one in ||]. The basic idea is to 
show that there is effectively no interaction between localized particles and low momentum photons. 
To make this idea explicit we write our Hamiltonian in a gauge different from the usual Coulomb 
gauge. 

To be precise, define 

A(x) = A(x) - A(0), (47) 
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which is well defined owing to the ultraviolet cutoff. The unitary operator that accomplishes this 
is U = exp 

b'EjliV^r ^(o)]. This is an 'operator- valued gauge transformation'. It commutes 
with A(x) for all x, but not with ojfc) or with Hf. 
Define 

b x (k, X) = Ua x (k)U* = a x (k) - iw x (k, X), (48) 
with w\(k,X) = xa(^)|^|~ 1 ^ 2 ^a(^) • Ylj=i x j- The transformed Hamiltonian H m is 

H m = UH^U* = ^ (p 3 + VZA(xjj) + |v^^ • B(xj) \+V + H f (m) (49) 

= E f "m(k)b x (k,X) b x (k,X) . (50) 

A=l,2 ^ K 

To estimate ||a A (fc)$ m ||, write 

a x {k)<5> m = U*a x (k)$ m - iw x (k, X)<5> m (51) 
where $ m = U<& m , and note that 

\\w x (k,X)$ m \\ < v^^|||X|* m ||. (52) 
It remains to estimate ||a A (/c)<3? m ||. By the Schrodinger equation for <3? m 



H m -E v (m,N) \ a x {k)$ m = H m ,a x {k) 



AT 



j'=i 

g kAe x (k) „ - ik . Xj % 

-Lu m (k)b x (k,X)$ m . (53) 

While this equation is correct, the derivation is somewhat formal. This is rigorously justified 
in Appendix [B]. 

Now add uj m {k)a x {k)<b m on both sides. Since E v (m, N) is the ground state energy and co m (k) > 
the operator H m — E v (m, N) + oJ m {k) has a bounded inverse R(u) m (k)) and hence 



a A (fc)<K m =2^R{u m {k))^te x {k) ■ + V^A{x 3 ))(l - e~ lk 



v 



V 



3=1 

o r—k A £ X (k) (54) 



- R(oj m (k))uj m (k)iw x {k, X)$ m . 
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For consistency, note that for \k\ > A w\(k,X) = and hence a x (k)& m = 0, i.e, for these modes 
<& TO is the vacuum as it should be for a minimizer. Since ||fl(o; m (A;))|| < uj m {k)~ l the norm of the 
last term is bounded by 



^f§lll*|5 m ||. 

To bound the norm of the first term we need to estimate 

N 

\\Y,R^m{k))e x {k) ■ ( Pj + y/ZA( Xj ))(l - e~ ikx ^ 



(55) 



N 






3=1 






N 


sup 


E( 


T?||<1 


3=1 




N 


sup 


£ 


77II < 1 


3=1 



^A( Xj ))R(u; m (k))r,f 



Next estimate the square of the first factor to get 



1/2 



JV 

End 

3=1 



e ik ' Xj )®m\\ 2 



1/2 



(56) 



ri,R(cj m (k)) 



N 



R(u> m (k))i] 



_i=i 

< a (77, R{u m {k))H^R{u m {k))ri) + 6 

< a (7/, R(uj m (k))rj) + (a^fm, iV) + 6) (77, R{u m {k)frj) 
,(A + 1) 



< C- 



\k\< 



for I A; < A 



(57) 



where a and b are independent of m. Since sup m<1 E v (m, N) < 00 the constant C is also indepen- 



dent of m. Finally the second factor in (56) is bounded by |/c||||X|$ m ||. The term containing the 
Pauli matrices in (54) is estimated similarly. In conjuction with ((54]) , (|5^)(^q) and ( |57|) this shows 
that 



\a x (k)$ m \\ < C\/a(A + l 



,1/2 



XA(fe) 



\x\$ r 



This, together with Lemma 6.2, proves the theorem. 



(58) 
□ 



Next we differentiate (|54| ) with respect to k. There is a slight problem with this calculation since 
the polarization vectors cannot be defined in a smooth fashion globally. We make the following 
choice for the polarization vectors. 

(k 2 ,-h,0) 



ei(k) 



(59) 



and 



e 2 (k) 



Aei(fc) . 



(60) 
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6.3. THEOREM (Photon derivative bound). Assume that there is binding, i.e., T, v (m, N) — 
E v [m,N) > 0. Assume that <I> m is a normalized ground state for the many-body Hamiltonian H^, 
m > 0. Then for \k\ < A and (fcj, k 2 ) ^ (0, 0) 

nv».»(t)M < ^Fra- < 61 > 

where Q is a finite constant independent of^> m , g, a, A, and depends on m only through the binding 
energy T, v (m,N) - E v (m,N) > . 



Proof. We differentiate 54 with respect to k and obtain 
V k (a x (k))$ m = 



2 v / a^K 1 (/c)) 2 — e A (fc) • + v / a^l(x j )) i ~^ 1/2 Jj $ m + 

2 v / a^K l (/c))V fc (£ A (A : )) • ^(p,- + ^A{ Xj )) { ~* Jj $ m + 



2 v / a^i?(o; m (A:))e A (A;) • + yfeA(xj))Vk f T^ 1/2 J $n 

j=i 



-V*. fi?(a; m (A;))a; m (A;)w A (A;,X)$ r 



The norms of the first and third terms can estimated precisely the same way as in (56) and 
([57]), and yields a bound of the form 

^ii(i+m)«U • (62) 

For the second term, a straightforward calculation shows that 

const 

|V fe £i(fc)| < — ===== for i = 1,2. (63) 
V + 2 

The last term is dealt with in a similar fashion as the previous ones. Using the steps in (^) and 
([57]), this leads to the bound 

l|V t (a iW )8 m ||< wl/2 ^ + ^ ||( 1 + |X|)8 m ||. ,64) 

The fourth term can be estimated in the same fashion to yield a similar result. 

Differentiating (|5l|) leads to the same estimate with 3> m replaced by <& m . This, together with 
Lemma 6.2, proves the theorem. □ 



As for the proof of Theorem 3J. , our somewhat formal calculations above are rigorously justified 
in Appendix pL 
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7 PROOF OF THEOREM |2H] 

The proof will be done in two steps. 



Proof. Step 1. The Hamiltonian has a normalized ground state $ m , by Theorem 4.1. Pick 
a sequence mi > m 2 > ... tending to zero and denote the corresponding eigenvectors by This 
sequence is a minimizing sequence for i?Q by Theorem |5,3| . Since is bounded there is a 

subsequence (call it again <F,) which has a weak limit Since H% - E v (0, N) > and by the 
lower semi-continuity of non-negative quadratic forms (in our case, Hq — E v (0,N)) 

0< ($,(HX -E v (0,N))$) < liminf($ 3 -,(^,5 / - E v (0,N))<f>j) = . 

Hence $ will be a (normalized) ground state if we show that ||$|| = 1 (i.e. &j —* $ strongly). 
It is important to note, however that if we write f>, = {$j 0, $4,1, &jn> ■••}; where <fr, n is the 
n-photon component of $j then it suffices to prove the L 2 norm-convergence of each &j in - The 
reason is the uniform bound on the total average photon number; see the remark after Theorem 



6.1| which implies 

^2 ll*J>f < const iV -1 . 



n>N 

-2 



Likewise, it suffices to prove the strong L convergence in the bounded domain in which \X\ < R 



for each finite R. The reason for this is the exponential decay given in Lemma 6.2, which is uniform 



by Lemma 5A. Finally, by Theorem |6.1| <&j^ n (X, fa, . . . ,k n ) vanishes if \fa\ > A for some i. So it 
suffices to show L 2 convergence for $j, n restricted to 

Q = {{X,fa,... ,k n ) : \X\ < R; \fa\ < A, i = 1, . . . , n} C R 3 ( Ar+?1 ) 

for each R > 0. 

Step 2. For each p < 2 and R > we show that $j n restricted to f2 is a bounded sequence in 
W llP {n). The key to this bound is fH]) and 

(a x (k)$j) n -i(X, hi, . . . ,k n -i) = T/n$j >n (X,k, fa,... , fc„_i) (65) 

where the arguments A, Ai, . . . , A n _i and the spin indices have been supressed. By the symmetry 
of %,n> (|65D, Holder's inequality and ( |6~i"| ) 



f dX f dfa... dk n £ \V ki $ jtn (X, fa, . . . , kn)\ P 

JB R J\ki\,...,\k n \<A i=1 

= n 1 ~P/ 2 [ dX [ dfa...dk n \V kl (a x {fa)$ j ) n - 1 (X,k 2 ,... ,k n )\ p 

JB R J\k 1 \,...,\k n \<K 

<C f dfal [ dX [ dk 2 ...dk n \V kl (a x (fa)$ j ) n -i(X,k 2 ,... ,k n f 

J\ki\<A \J B R J\k 2 \,...,\k n \<A 

<C I dfaWV^a^fa^jW 13 < const (66) 

J\ki\<A 



P/2 
|2 1 
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independent of j. The constant C depends on all the parameters, but is finite because \k%\ < A in 
the integration. Similarly, by Holder's inequality 



\ X (\X\ < R)V x *j,n\\$ < C\\ X (\X\ < R)V X ^ 



n\\2 



N 



< C(^,5>^. )P /2 
i=l 

which is uniformly bounded by Lemma ([A. 5 ). 

Since the classical derivative of a x (k)& m is not defined in all of SI one has to check that the 
weak derivative coincides with the classical derivative a.e.. Because of our definitions (|59| ) and (|60|), 
the classical derivative is not defined along the 3-axis. 

One has to show that 

/ di^ j)n = lim / ditp^ j!n = - lim / V<%$j> i = I, ■ ■ ■ , 3(iV + n). 
Jn £ -^°Jn e £ ^°Jn E 

for any test function ip £ C£°($l). Here Q s is SI with an e cylinder around the 3-axis removed in 

each fc-ball. The first equality is trivial; it is the second equality that has to be checked. This 

amounts to showing that the boundary term, coming from the integration by parts vanishes in the 



limit as e tends to zero. But this follows immediately from Theorem 6.1 



This shows that $j >n as a function of all its 3(N + n) variables, is in the Sobolev space W 1,p (£l) 
and that sup^ ||3>j,n||wi.i>(n) < 00 • Since $ n j converges weakly in L 2 (Q) it converges weakly in 
L P (S7) and since the sequence is bounded in W 1,P {Q), V$ n j converges weakly to ^7*& n . 

The Rellich-Kondrachov theorem (see |22j Theorem 8.9) states that such a sequence converges 

strongly in L 9 (S1) if 1 < q < [3p(N + n)/3(N + n) — p\. The boundedness of f2 is crucial here. For 

our purposes we need q = 2, and hence we have to pick p such that 

2-3(iV + n) . . 

2>?> 2 + 3(iV + n) < 67) 

which is possible for each N and n. We conclude that <&j, n — * &n strongly in L 2 (S7) as j — > oo, for 
each n and R. This proves the theorem. □ 

Remark: Theorem |6.3| essentially says that the derivative is almost, but not quite in I? . For 
high dimensions, 3(iV + n), the required p is as close as we please to 2 if we require q = 2, but p = 2 
is not allowed. The way out of the difficulty was to prove a uniform bound on the number operator 
and use this to say that that it suffices to prove strong convergence for each n separately. With 
n then fixed, it is possible to find a p < 2 that yields q = 2. It is, therefore, crucial to have the 
derivative in every LP space with p < 2. The resolution of the problem of the infrared singularity 
is thus seen to be a delicate matter. 



A Appendix: LOCALIZATION ESTIMATES 

In this appendix we collect a few facts which we use several times in this paper. Generally we worry 
about localizations of Hamiltonians in configuration space. While this is standard for Schrodinger 
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operators it is somewhat more complicated in the presence of the radiation field. This is chiefly 
due to the problem of localization of photons. 

We begin by stating a few well known facts about partitions of unity. Let j3 denote one of the 
2^ subsets of the set of integers 1,2,..., N. Its complement is denoted by (3 C . As shown in |2C] 
there exists a family of smooth functions jg having the following four properties. 

(i) 

13 

(ii) For (3 ^ {1, . . . , N} (including the empty set) the jg's are homogeneous of degree and live 
outside the ball of radius R centered at the origin and 

supp jg C {X : min (\xi — Xj\, \xj\) > c\X\} , (69) 

where C is some positive constant. 

(iii) In the case where (3 = {1, . . . , N}, jp is compactly supported. 
Corresponding to these electron localizations we define photon localizations. 
For given (3 ^ {1, . . . , N} consider the function 

gi (y;P,X)=U je0 c fi- x (1LZEi)\ (70) 

where x 1S a smoothed characteristic function of the unit ball. Define 52(2/; P, X) = 1 — g\(y, (3, X). 
In the variable y, the function g\ is supported away from the particles in (3 C while g 2 lives close to 
the particles in (3 C . Next define, for i = 1,2, 

ji {y;(3,X)= - 9 MM1 (71) 
Vgi(y,f3,Xy+g 2 (y;[3,Xy 

Certainly j\ + j\ = 1 and a simple computation shows that 

,_ . . const. 

|Vji| < — — . (72) 

In the case where [3 = { 1 , . . . , N} the construction of ji and j'2 is similar to the above one except 
that the function g\ depends on y, is equal to one in a neighborhood of the origin and is compactly 
supported. 

With the help of j\ and j% the photons can now be localized as was done in Section 4. Let 
Ug(X)J- — > J- <S> T be the corresponding isometric transformation, i.e., the one that is defined via 
the relation 

Ug{X)a*{h)U* 3 {X) = a *{j l h)® l + l®a*{j 2 h) . (73) 

The tensor product indicated is a tensor product between Fock spaces. 

We denote by Hg the Hamiltonian of the form (Q) with photon mass, but only for the particles in 
the set (3. More precisely this operator acts on L 2 (M 3 '^') (g> T . By H@ c we denote the Hamiltonian of 
the form (Q) with photon mass, but only for the particles in the set (3 C where the interaction with the 
nuclei has been dropped. This operator acts on L 2 (R 3 '^ C ') T In particular we keep the interaction 
among those particles. In the case where j3 = {1, . . . , N} the Hamiltonian H@ c = Hf{m). 
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A.l. LEMMA (Localization of Hamiltonian) . For every (5 

jpHjp = U* p (X)jp [h p ®1 + 1® jfsUfs(X) + . (74) 

For (3 / {1, N}, o(l) -> as first R -> oo and then P -> oo. If (3 = {1, ... , N} then o(l) -► 
as P ^ oo for every fixed R > 0. 

Proof. Our immediate aim is to compare the field energy iij with the localized field energy 
Up(X)[Hf ® 1 + 1 ® Hf]Up(X). For simplicity the various indices are supressed and Up(X) is 
replaced by Up. The variable X plays no role here. Pick an orthonormal basis {gj}JL l of L 2 (IR 3 ) 
in H 1 ! 2 ^). States of the form 

C = const. a\ n (g h )- • • (g ik )\0 > (75) 

where A; is finite, form an orthonormal basis in the Fock space. The field energy acts on such states 
as 

k 

H fC = E a K ^n) ■ ■ ■ a K 3 i^ij) ■ ■ ■ a* Xik ( 9lk )\0 > . (76) 

Thus, we have that 

k 

H fC = U P E c k ^n)' ■ ■ 4, (^) • • • 4 fc (9i k )Up\0 > (77) 

3=1 

and 

H f = U*p[H f ®l + l®H f ]Up + E f (78) 

where the error Ef is given by 

k 

^ 1 3 3 (79) 

■■■< ifc (^)^|o> . 

Thus Ef is given by the operator (the A's are omitted) 

E f =Up^2[a*(\jiM9k) ®1 ' + 1 '® a*(\j 2 ,Lu}g k )}[a {j igk )®l + l®a (j 2 gk)]Up . (80) 
k 

The expression for the operator Ef does not look hermitian but it is, remembering that Up is a 
partial isometry. Standard estimates lead to 

| (*,£,*) | < (H^lll + \\[u,j 2 }\\) (¥, [AA+ 1] tt) . (81) 



where is the number operator. 
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Here ||[ji,cj]|| denotes the operator norm associated with the kernel This norm can be 

estimated using the formula 



= \ji,u 2 )- +w 2 [ji, -] . 



(82) 



Recalling the definition of ji, the operator norm of the first term is easily seen to be bounded by 
a const./ P. Likewise, the second term, using the formula 

1 1 f°° 1 dt 

(83) 



yV + m 2 n Jo t + P 2 + m 2 ' 

can be estimated by const./ P. The term ||[j2,u;]|| is estimated in a similar fashion. The estimate 
(|3l|) immediately shows that for a general state $ we have that 



const. 



|($, [H f -U$[H f ®I + I®H f ]Up] $)| < — — ($,AA$) . 



(84) 



Since the photons have a mass we can estimate the number operator in terms of the field energy. 
The field energy is relatively bounded with respect to the Hamiltonian, i.e., Hf < aH^ + b for 
some positive constants a and b, and thus we obtain 



const. 



|($, [Hf -U;\Hf®l + l® H f ] Up) n < -p^(*, [aHl + b\<S> 



(85) 



Note that this estimate had nothing to do with the electron, in particular the x-space cutoff is not 
present in the calculation. 

Next we have to compare YljLiiPj + \/a A(xj)) 2 with 



^2(pi + ^/aA(xi)) 2 + ^(Pj + ^faA( Xj )f 

ie/3 je/3 c 



This time the X-space cutoff is important. We would like to estimate the difference 
J^[(K + V^A( Xi )) 2 - U;{X){ Pi + V^A( Xi )) 2 <g> TUp(X)] j p 



i€P 



+ jp [(p* + V^A( Xl )) 2 - u;(x)i®( Pi + ^A{ Xl )) 2 u p {x)\ j p . 

ief3 c 

It suffices to treat the first term, the other is similar. It can be easily expressed as 



(86) 



Jp 



^2(pi + VaA( Xi ))Qi + Q l ( Pi + y/aA(xi)) - Q\ 

i&P 



(87) 



where 



Qi = Pi + V^A{ Xl ) - U*JX)(pi +^A(xi))®lUp(X) . 



Using the form boundedness of the kinetic energy with respect to the full Hamiltonian, we have 



N 



(89) 
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for positive constants a and b.. Thus, using Schwarz' inequality it suffices to show that 



11^*11=0(1) for iefi, 
as R (the localization radius for the electrons) tends to infinity. Denote by 

= (27T)- 3 / 2 

Explicitly, Qi is given by 

Vi -U* p (X) Pi ®lUp{X) 
+ Up{X) 



\k\<A u{k) 



(90) 



(91) 



+ U%{X) 



J2 a x(.[h ~ 1]^*) + a x {j 2 h 

. A 



Up{X) 
Up(X) 



(92) 



and it suffices to estimate each of these terms separately. Each of the last two terms can be brought 
into the form 



\\U* p {X)a*{f)®lUp{X)j^\\ 

where / is one of the functions 

[j 1 (y,f3,X)-l]hl Xj (y) or j 2 (y, (3, X)h^. (y) 
The terms (|93|) are estimated by 



j€0 



X J 



respectively 



A" 



(93) 



(94) 



(95) 



(96) 



In both formulas the index j is in (3. The function jp lives in the region where \xi — Xj\ > cR for 
i € (3 and j G (3 C . The function j% — 1 (and likewise j 2 ) is not zero only if \y — xj\ < P for some 
j £ (3 C . Thus, jp(X)(Ji — l)(y) and jp(X)j 2 (y) are nonzero only if \y — Xi\ > cR — P. As cR — P 
gets large only the tail of the function h x contributes to the integral which can be made as small 
as we please. The number operator is bounded by the field energy times 1/m which in turn is 
bounded by the full energy. 



To estimate the first term in (92) we calculate 



Pi -U%(X) Pi ®lUp(X) = 
U;(X)Y,[a*{[PiJi}9k)^l + l0a*([ Pl ,j 2 }g k )} 



(97) 



[a (jig k ) ®T + T®a (j 2 g k )j Up(X) . 
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Note that the tensor product in the first line is different from the second. In the first the identity 
acts on L 2 (M 3 I^ C I) ® T ® T while in the second ® indicates the tensor product of the Fock spaces 
only. The functions g k indicates a basis of L 2 (M 3 ). The operators U^(X) and Up(X) have unit 
norm. Thus 

\\Up(X) [a*(\Pi,3l}9k)®I + l®a*(\pi,j 2 }g k )] x 

k 

[a (hg k )®l + l®a (j 2 g k )} Up(X)*\\ (98) 
<(ll[Pi,Ji]|| + ll[Pi,j2]||)||v<Ar+T*|| , 

where || • || indicates that the operator norm has been taken. The norms of the commutators are 
of the order 1/P and hence vanish as P — » oo. Since the photons have a mass we can estimate the 
number operator in terms of the field energy. 

Similar consideration apply to the f3 c term in (86). The only difference is that instead of ( |95| ) 
and (^6|) we have 



respectively 



su V {jfs{X)\\hh x x .|| 2 }V(*, (Af + (99) 
x 3 



sup{j>PO||L72 - l]h x x . \\ 2 }^f(^W+W) , (100) 

x 3 



with j £ p c . Again this terms tend to zero as P — » cxd. The proof for the case where (3 = {1, . . . , N} 
is similar but simpler since the operator Up does not depend on X. 

Finally, we have to compare the a ■ B term with its localized counterparts. The estimates 
are similar to, but much easier than the estimates for (p + yfa A(x)) 2 and are omitted for the 
convenience of the reader and authors who, by now, are exhausted. 

□ 



A simple consequence of Lemma A.l is the following. 

A.2. COROLLARY. Let cp be a smooth function on R 3N such that jptj) = for /3 = {1, . . . ,N}. 
Thus, (j) depends on R. Then, as operators, 

cpHcf) > (S y (m, N) + o(l)) <\> 2 . (101) 

Here, E v (m,N) = mm 1 < N/<N (E v (N') + E°(N - N')) and o(l) vanishes as R -> oo. 

Proof. By the IMS localization formula we have that 

<j>H4> = <PjpHjp<t> -4> 2 Y I v ^l 2 . ( 102 ) 

p p 

where the second term goes to as R — * oo. With our assumption on cp only the sets (3 with (3 C ^ 
contribute. From Lemma |A.1| we get that 



<PH4> = Y J U*p(X)^jp [Hp®T + T®H pc \jp<l ) Up{X) + o{l) (103) 
& 
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as first R — > oo then P — ► oo. Certainly .ffg > E v (m, \j3\) and > E°(m, \(3 C \) from which the 
statement immediately follows. □ 

A. 3. LEMMA. Let ^> n be a normalized sequence in 7i whose energy is uniformly bounded and 
such that for any with finite energy, 

*) -+ , and H$) . (104) 

T/ten 

(fe*n, C/*T ® iW0A) - . (105) 
.Here J7 is i/ie FocA; space localization Up that corresponds to (3 = {1, ...,iV}. 
Proof. Since the energy of *S> n is uniformly bounded we also know that 

(¥ n ,-ff°(m)tf„)<C7. (106) 

is uniformly bounded 

Let us describe the operator I ® P2U in more detail. Recall that 

t/a*^) • • • a*(h ik )\0 >= c*(h h ) ■ ■ ■ c*(h ik )U\0 > (107) 

where |0 > denotes the vacuum vector in Fock space and U\0 >= |0 > &)|0 >. Hence, using the 
definition of c*(h), we find that 

J ® P 2 Ua*(h h ) ■ ■ ■ a*(h ik )\0 >= a*(h h h) ■ ■ ■ a*{h ik ji)\0 > ®|0 > . (108) 

The projection P2 annihilates the photons in the second factor. In other words, the operator 
I ® P2U when acting on a state 

^ = {M/°,* 1 (y 1 ),M' 2 (y 1 ,y 2 ),---} (109) 
produces the localized state F(ji)^ (g) |0 > where 

r(i 1 )^ = {^°,i 1 (y 1 )vi/ 1 (yi),ji(yi)ji(y2)^ 2 (yi,2/2),---} ■ (no) 

It follows that 

= ||rO'l)0B*n|| 2 

Next, we show that ( |106| ) implies that 

r(il)0i?^n -> . (Ill) 
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To achieve that we note first that on account of the positive mass we have that (\E r n ,A/"\E , n ) is 
uniformly bounded. Since ^ n is of the form 

{* n ,* 1 n (X,y 1 ),y 2 n (X,y 1 ,y 2 ),---} 

we know that Ylk>M (^n> ^n) — const/M. It is therefore sufficient to prove ( [LllD for each function 

*%(X, yi ,--- ,y M ) ■ 

From the lemma below we learn that 

^2^ n ,pp n ) (112) 

3=1 



is uniformly bounded. Thus, we can write ( |111|) as 

r(j 1 )^(i + + -ff/) _1 ^ 2 (i + J>? + ^/) 1/2 *« • ( 113 ) 

i=i j=i 
which vanishes as n — > oo since ||(1 + X^=iPj + -^V) 1 ^ 2 ^™!! * s uniformly bounded and since 

r(ji)<Mi + 5> 2 + tf/r 1/2 

is compact on every finite particle subspace. Compactness follows from the fact that for continuous 
functions / and g vanishing at infinity the operator f(iS7)g(x) is compact. □ 

A. 4. LEMMA (Bound on A(x) 2 ). For each x G M 3 and ultraviolet cutoff A write A(x) = 
D{x) + D*(x) where D contains the annihilation operators in A(x) and D* the creation opera- 
tors. Similarly, write B{x) = E(x) + E*{x). As operator bounds 

H f > ——D*(x)D(x) 
1 ~ 8vrA v ; v ; 

Hf + T > ^(x)E*(x) 

Hf + — > -J—B(x) 2 . (114) 
1 16 ~ 32vrA 3 y ' y ' 

Proof. We write A(x) = D(x) + D*{x) with D(x) = £ A Li <A |A;| _1 / 2 £a(^) exp[iJfc ■ x]a x (k)d 3 k. 
There are thus four terms in A(x) 2 . Using the Schwarz inequality, the (DD) term can be bounded 
above by (D*D)/2 + (DD*)/2. On the other hand, (DD*) = (D*D) + T, where V is the commutator 
/ 2/\k\ = 4vrA 2 ; the factor 2 comes from the two polarizations A = 1,2. Altogether, we obtain 

A(x) 2 < 4D*(x)D(x) +4ttA 2 . 



MGELML/3/2/01 



28 



Finally, we use the Schwarz inequality again to obtain 

[Mk)a* x (k)d 3 k [ h x (k)a x (k)d 3 k<Y^ [ \hx(k)\ 2 /\k\d 3 k £ I \k\al(k)a x (k)d 3 k. 



A 



In our case, h\(k) = £\(k) exp[ik ■ x]/y/\k\, so J2\ / I^a(^)| 2 /|^M 3 ^ = 87rA. 

For B = curLA, replace T by 27rA 4 and replace |/ia(&)| by \/\k\. □ 



As a corollary of Lemma [A.4| we have the following. 

A. 5. LEMMA (Bound on (p + A(x)) 2 ). For any e > there are constants 5(e) > andC(e) < 
oo such that 

N N 

J2{(Pi + V^A( Xj )) 2 + 9 -^a 3 ■ B(xj)] + eH f > 5(e) J> 2 - C(e) . (115) 
The constants 5(e), C(e) depend on a, g, A, N. 

Proof. In addition to Lemma [A. 4 , use the facts that for any < \i,v < 1, (pj + y/a A(xj)) 2 > 
(1 - fi)p 2 + (1 - l/fj)aA(xj) 2 and 2a j ■ B( Xj ) > -vB(xj) 2 -\jv. □ 



B Appendix: VERIFICATION OF INFRARED BOUNDS 

The proofs of the infrared bounds in Section |6| are somewhat formal. In particular, we carried out 
the calculations tacitly assuming that a\(k)<& m (which is itself only defined for almost every k) is 
in the domain of the Operator H m . One can actually prove this when H m is self-adjointly realized 
in terms of the Friedrichs' extension and thereby make all the formal computations in Section ^ 
rigorous. Instead of doing so, we give here alternative proofs of the Theorems in Section [6| which 
avoid any reference to a domain of H m . All the arguments can be carried out on the level of 
quadratic forms. 

We recall that H m is the Hamiltonian after an "operator- valued gauge transformation". 
Our remarks here about quadratic forms in relation to H m could just as well be applied to 
itself. 

In order to keep the notation simple, we give the proof of the infrared bounds for the case of 
a single charged particle (N = 1) with no magnetic moment, i.e., g = 0. There is no difficulty in 
deriving these bounds for the general case. 

Denote by S the set of all finite linear combinations of vectors that are products of C^°(1R 3 )- 
functions and states in T that have only a finite number of photons. It is well known that this set 
is dense in Ti, and that the quadratic form ('I', \&)_|_ := (\&, (H m — E v (m, 1) + 1)^) is defined for all 
^ in S and is bounded below by H^H 2 . Hence this quadratic form is closable and the closure of S in 
this inner product is a Hilbert space Q(H m ) with inner product (•, -) + and norm \\^\\+ = \J (\E', ^)+- 

An eigenfunction <3? m of H m in the weak sense is a vector in Q(H m ) such that 



(116) 
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for some real number e and for all \& S Q(H m ). It is in this sense that we proved in Section [| that 
a ground state exists for the model with massive photons. (This implies that $ m is in an eigenstate 
of the Friedrichs' extension of H m ). 
Define the smeared operators 



<*(/) = F / a x (k)f(k,X)dk , (117) 



F 

A 

-2/ T n>3. f p2 



where f(k,X) is any function in L (R ;C ). It is not difficult to show that a(/)$ m is in the form 
domain of H m . To this end define a R (f) = R[Af + i?] _1 a(/). Here R is some large real number 
(which we eventually take towards infinity) and J\f is the number operator. 

It is straightforward to see that a R (f) and a R (f) are bounded operators on Q{H m ) for every 
R > 0, i.e., 

\\a R (f)n+ < c(R)\\n+ , (us) 

and similarly for a* R (f). 

Generally, the constant C(R) tends to oo as R tends to oo. For an eigenfunction of H m , however, 
this is not the case. Simple but tedious commutator estimates reveal that for any eigenfunction 
$ m there exists a constant C independent of R such that 

(a B (/)* m , a R (/)$ m )+ < C(a R (f)$ m , o B (/)* m ) . (119) 

The point is that (or(/)$ to , a R (f)$ m ) + = (a* R (f)a R (f)$ m ,$ m ) + plus terms that are uniformly 
bounded in R. By the previous statement we know that a R (f)a R (f)$> m is in Q(H m ) and hence 

(a R (f)a R (f)$ m , $ m )+ = e(a* R (f)a R (f)* m , £ m ) = e(o fl (/)* m> o B (/)* m ) . (120) 

The last expression, however, is bounded uniformly in R, since the condition $ m E Q(H m ) implies 
that the expectation value of the field energy in <I> m is finite which in turn bounds the last expression 
in ( |12C| ). Here we use the fact that the photons have a mass. 

From this it follows easily that for a subsequence of i?'s tending to infinity, aji(f)^ m has a weak 
limit in Q(H m ). Since a R (f)$> m — > a(/)$ m strongly this shows that a(/)$ m G Q(H m ). 



Proof of Theorem EO . We shall use the abreviation 

^ f ...dk = J ...dk . (121) 

A ^ 

For our special choice of gauge 

= o(G*) + a*(G*), i = 1, 2, 3 , (122) 

where we set 

G\(k,x) =e\(k)\k\- 1/2 (e ik - x -l) X A(k) . (123) 
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Next, pick any ^ in 5 and calculate (recalling the definition of w in Section y equation fl4q)) 
(¥, (H m - E v (m, l))a(/)$ m ) = -2(¥, (/, G)(p + A)* m ) - (tf , a(w/)* ro ) + (/, ^)5 m ) 



(124) 



with cj(A;) = \/k 2 + m 2 . This extends, using an approximation argument, to all ^ £ Q(H m ) and, 



in particular, to a(/)<3? m . Here we note that, on account of Lemma A. 4 and the assumption on the 
potential, \I> £ Q(H m ) implies that (p + -A) 1 ^ £ H. Hence 

< (a(f)$ m ,(H m -E v (m,l))a(f)$ m ) (125) 
= -2(a(/)$ m , (/, G)(p + - (a(/)* mj a(w/)$ m ) + i(a(/)* m , (/, c^)$ m ) , 



which yields the inequality 

(a(/)$ m , a(w/)* m ) < -2(a(/)$ m , (/, G)(p + I)$ m ) + i(a(/)* m , (/, uH$ r 



(126) 



for all / in L 2 (R 3 ;C 2 ). Pick / of the form uj(k)~ 1 / 2 q(k, X)gi{k, A) where <7i is an orthonormal basis 
of L 2 (1R 3 ; C 2 ) and q(k, A) a bounded function. Summing over this basis, we get on the left side of 

S 

^(a(^ 1/2 ^)$ m ,a(u; 1/2 ^5 m ) = ^ \q(k, \)\ 2 \\a x (k)$ m \\ 2 dk , (127) 

i 

and on the right side 

-2(a(a;- 1 |( ? | 2 G)i m , (p + A)l m ) + i(a(|(?| 2 u;)$ m , $ m ) . (128) 

Hence 

\q(k, A)| 2 ||a A (£0$ m || 2 d£; < -2(a(u;- 1 k| 2 G)$ m , (p + I)$ m ) + i(a(|g| 2 ^)$ m , $ m ) . (129) 

The right side can be written as 
\q(k,X)\ 2 



{a x (k)<5> m ,G x (k)(p + A)$ m )dk + ill \q(k, \)\ 2 (a x (k)<S> m , w x <S> m )dk , (130) 



and, applying Schwarz's inequality, this is bounded above by 

1/2 



(k,X)\ 2 \\a x (k)^ m \\ 2 dk 

Lo(k)- 2 \q(k,\)\ 2 \\G x (p + A)$ m \\ 2 dk 



(131) 



'|fe|<A 

Hence we obtain the bound 

|2l|„/UAA 112. 



1/2 



+ 



|g(fc,A)| 2 ||«; A $ m || 2 dfe 



fe|<A 



1/2- 



|g(M)rlKfc)*mlrdfc < 



(132) 



w(A;)- 2 | 9 (A;, A)| 2 ||G A (A;) • (p + A)$ m \\ 2 dk + 2 | g (A;, \)\ 2 \\w x <S> m \\ 2 dk 

\k\<A J\k\<A 
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Since, div^ G\ = we have that G\ ■ (p + A)<& m = (p + A) ■ G\<& m . Moreover, (p + A) 2 is relatively 
form bounded with respect to H m . But, as in the proof of exponential decay (Lemma |6.2[) , we have 



for each i = 1,2,3 



(G\$ m , (H m - E v (m, l))G\*m) = ($m, I V x Gl| 2 $ m ) 



and we arrive at the bound 



/(£;,A)| 2 ||a A (A0<i> m || 2 d£;<C 



fe|<A ^(k) 2 



IGa^H 2 + \\\V x G X \<5> m \\ 2 + Uj{kf\\w X <$>r 



(133) 

dA; , 
(134) 



where C is some constant independent of m. Since q{k, A) is arbitrary we obtain for almost every 
k and each A that 



a x (k)$ m \r < Cu;(k)- 2 \\G x $ m \\ 2 + \\\V x G x \$ m \\ 2 + u J (ky\\w x <5> m \\ 2 X A(k) 



The right side is bounded by 



C ~ 

rn\\\ x \^m\\ 2 XA(k) , 



which is finite on account of the exponential decay of <3? m . 



Proof of Theorem \6.3[ . First some notation: For any function f(k) define 

(A h f)(k) = f(k + h)-f(k), 



and 



A_ fc o(/) = a(A h f) 



Returning to ( |126[ ) with / replaced by Ahf we have 



(A_ /l o(/)f m ,a(wA A /)$ m ) < 

2(A_ A a(/)* ro , (A h f, G){p + A)$ m ) + i(A_ fc a(/)* ro , (A h f, ujw)$ m ) 



which can be rewritten as 



(A_ h a(/)$ TO , A- h a(wf)$ m ) 
< (A_ fe a(/)$ m , o((A fc w)/(- + fc))* m ) 

- 2(A_ ft a(/)$ m , (/, A_ h G)(p + i)$ m ) + i(A_ fc a(/)* ro , (/, A- h (uw))$ r 



(135) 

(136) 
□ 



(137) 



(138) 



(139) 



(140) 



Notice that without the first term on the right of the inequality sign, the structure of this inequality 
is the same as (|126| ), except that, of course, A_/ l a(/) plays the role of a(f), A_^G plays the role 



MGELML/3/2/01 



32 



of G and A_j l {ujw) plays the role of cow . Thus, without this term we would obtain immediately 
the estimate analogous to (|134j) , 



< C 



(fc,A)| 2 ||(A_ h o A )(fc)* m [| 2 dfc 

k(M)l 2 



u(kf 



{A^ h a x (k)<$> m ,a x (k - h)$ T 
which, by Schwarz's inequality, is bounded above by 

1/2 



A_ h G x $> m \\ 2 + \\\VxA_ h G x \$ m \\ 2 + co 2 \\A_ h (ujwx)<S> m \r M . (141) 



(142) 



The remaining term in equation ( |14C| ), after summing over the functions qgi/y/uJ, turns into 

k(M)l 2 



u(k) 



-(A h oj)(k - h)dk 



\q(k,\)\ z \\A_ h a x (k)<S> m \\ 2 dk 



a x {k - h)^ m \\ 2 lq{ )^f \A h u;(k ~ h)\ 2 dk 



1/2 



• (143) 



This, together with (141), yields 

\q(k,\)\ 2 \\(A_ h a x (k)$ m \\ 2 dk 
l(k,X)\ 2 ^ 



< C 



u(kf 



A_ fc G A $ m || 2 + \\\V x A^ h G x \$ m \\ 2 + u]{k) 2 \\A_ h {uw x )$ m \\ 2 
q(k,X)\ 2 



dk 



+ C 



L0{kf 

Again, since q is arbitrary we obtain for every fixed A 



\a x (k - h)$ m \\ 2 \A h v(k - h)\ 2 dk . (144) 



{A_ h a x {k)$ m \\ 2 < 



2 - C 



to{k) 2 L 



IA^^Ga^H 2 + IHV^A^GAl^mll 2 

+ u(k) 2 \\A_ h (uw x )^ m \\ 2 + \\a x (k-h)^ m \\ 2 \A h u(k-h)\ 2 ] . (145) 



Combining this with (|136j ) we get 



\\(A_ h a x (k)<5> m \\ 2 < 



2 - C 



u>{kf 



\A_ h G x $> m \\ 2 + \\\V x A„ h G x \<$> m \\ 2 + L0{k) 2 \\A_ h (L0W X )$>r 



+ 



c 



\x\$ m \\ 2 X A(k - h)\A h u(k - h)\ 2 . (146) 



u(k) 2 \k-h\ 

The polarization vectors defined in (|59|), (|60|) , are differentiable away from the 3-axis. The same 
straightforward estimates as in Section [6| lead to 

1 1 



\\(A_ h a x (k)$ m \\ 2 < C 



\h\ 2 \\\x\<S> m \\ 2 (147) 



\k\{k 2 + k 2 2 ) \k-h\((k 1 -h 1 ) 2 + (k 2 -h 2 ) 2 )_ 

which hold for all \k\ < A and small \h\ with a constant C that is independent of m. 

Next, we observe that for k ^ fixed, there exist a sequence of h values, say hi, tending to zero 
so that h (A-h e a x {k)<& m converges weakly to some element Vj(k) which satisfies the estimate 
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Here ej is the j-th canonical basis vector. Next we identify Vj(k) as the weak derivative of a\(k)& m . 
This weak derivative, by definition, can be computed via the expression 

-(*,a(dj<l>)$ m ) (149) 

where *Z> is any state in H and (f> is any test function in C^°(IR 3 ). Clearly the above expression 
equals 

\\xn-j(^^. hej a x {k)^ m )4>{k)dk . (150) 

But along the sequence hi 

Um^(^,A_ hie .a x (k)^ m )cf>{k)dk = f ^, Vj (k))(f>{k)dk , (151) 

which identifies Vj(k) as the (negative) weak derivative of a\(k)Q m . □ 
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